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Abstract. Open-shutter RHESSI observations of 3-15 keV X-rays are found to 
exhibit active region transient brightenings and microflares at a rate of at least 
10 per hour occurring even during the periods of lowest solar activity so far in the 
mission. A thermal component fitted by temperatures of 6-14 MK dominates from 3 
keV to about 9 keV, but can be traced up to 14 keV in some cases, and has an average 
duration of 131(±103) seconds at 7-8 keV. The duration increases with decreasing 
photon energy. The peak count rate defined by cross-correlation is delayed at low 
energies. The temperature peaks early in the event and then decreases, whereas the 
emission measure increases throughout the event. The properties are consistent with 
thermal conduction dominating the evolution. In some of the bigger events, a second 
component was found in the 11-14 keV range extending down to 8 keV in some cases. 
The duration is typically 3 times shorter and ends near the peak time of the thermal 
component consistent with the Neupert efi'ect of regular fiares. Therefore the second 
component is suggested to be of non-thermal origin, presumably causing the beam- 
driven evaporation of the first component. The two components can be separated 
and analyzed in detail for the first time. Low-keV measurements allow a reliable 
estimate of the energy input by microflares necessary to assess their relevance for 
coronal heating. 



1. Introduction 



Plasma of many million degrees and unknown origin has been reported 
from solar active regions (e.g. Brosius et al, 1996). Observers of the full 
Sun find a power-law distribution for the differential emission measure 
(DEM) decreasing with temperature. The higher the level of activity, 
the larger is the coronal emission measure and the flatter its temper- 
ature distribution. In the presence of large active regions, it has been 
observed to extend to temperatures higher that 10 MK (e.g. Peres et al., 
2000, and references therein), but is not well-observed beyond 10 MK. 
On top of the power-law distribution, flares produce transient peaks. 
The DEM of flares has a peak at 6(ibl) MK for A2 flares increasing up 
to 23(±3) MK for XI flares (Feldman et al, 1996). At any phase, the 
flare DEM is quite narrow and covers less than a decade in temperature 
(Reale, Peres and Orlando, 2001). From observations unresolved in time 
and space thus emerges the question whether the high temperature 
plasma apparent outside flares is just the remnant of previous larger 
flares, is heated by a stationary process or results from a superposition 
of unresolved small flares. 
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Small flares in the 3-8 keV range lasting a few minutes have first 
been observed by HXIS on SMM (Simnct et al, 1989). Shimizu (1995) 
found microflares (active region transient brightcnings) in Yohkoh/SXT 
observations sensitive to 4-7 MK. He reported events having energies 
as low as a few lO^^ergs and a correlation of peak flux with source 
size (loop length) and density. Watanabe et al. (1995) have observed 
the full Sun in the high-temperature Sulfur XV lines using the BCS 
instrument on board Yohkoh, indicating brightenings in excess of 10 
MK. Note in contrast that the thermal energy content of micro-events 
in the quiet corona is more than two orders of magnitude smaller, and 
temperatures of only 1.0-1.6 MK are reported (Krucker et al, 1997, 
and reviewed by Benz and Krucker, 2002). 

Lin et al. (1984 and references therein) have discovered small HXR 
flares produced by >20 kcV electrons. Comparing the >25 keV channel 
of LAD/CGRO with microflares observed in the 8-13 keV channel of 
SPEC/CGRO, Lin, Feffer and Schwartz (2001) reported many low- 
energy events that have no counterpart at high energies. If interpreted 
as emission of the same non-thermal electron population, it would sug- 
gest that the spectrum steepens or cuts off between 10 and 25 keV. 
Thus the nature of the low-energy events and their relation to the 
non-thermal events at high energy remain questionable. 

In this paper, we present microflare observations by the RHESSI 
satellite (Lin et al, 2002) at unprecedented spectral and temporal 
resolution of photons above 3 keV, highly sensitive to the highest 
temperatures reported for active regions as well as to the lowest ener- 
gies of non-thermal photons expected from microflares. First RHESSI 
results on microflare imaging are presented by Krucker et al. (2002). 
Here wc explore the spectral and temporal characteristics of microflares 
and concentrate on the question whether the emissions are thermal or 
non-thermal. 



2. Observations and Results 

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) 
was launched into orbit on 5 February 2002. The germanium detectors 
register photons in the energy range from 3 keV to 17 MeV with 1 
keV resolution at low energies (Smith et al, 2002). Imaging is achieved 
by nine absorbing grids modulating by satellite rotation. It allows to 
reconstruct the full Sun image with a resolution of 2" at low energies 
(Hurford et al, 2002). At low count rates, the spatial resolution is 
degraded (Saint-Hilaire and Benz 2002). RHESSI observes the Sun 
full time interrupted only by satellite night and the South Atlantic 



microflares.tex; 2/02/2008; 0:50; p. 2 



Microflcires in solar active regions 



3 



Anomaly. Energetic electrons sometimes cause spurious counts when 
the satellite is at high geographical latitudes. We use standard RHESSI 
software and calibrations updated to mid-August 2002. 

As the energy, time and detector information of each photon is 
registered, the large number of low-energy photons during large flares 
causes pile-up in the detector. To bring it down to a manageable level, 
one or two shutters are usually put in front of the detectors to reduce 
the count rate at low energies. Early in the mission, the shutters were 
open only in relatively quiet times, when no large flares were expected. 

Here we have selected RHESSI observing times when no shutter 
was absorbing the low-energy photons. During these times the full flux 
of low-energy photons is available. Thus these times yield maximum 
sensitivity for low-cncrgy studies. We have found some 50 open-shutter 
orbits in the months of March, April and May 2002. Uninterrupted 
intervals at low solar activity and low background were selected for 
further analysis. In particular, we required (i) a count rate less than 
600 cts per second and per detector in the 3-12 keV range and (ii) 
a constant background in the 100-300 keV range. A total of seven 
intervals satisfled these requirements. The GOES level in all of them is 
below C, and all flares are below GOES class A9 after subtracting the 
background. 

2.1. Light Curves 

Figure 1 displays RHESSI full Sun observations at low energies during 
one-orbit intervals. The low peak count rates justify neglecting the pile- 
up effect of photons in the detectors. The photons have been integrated 
in one keV bins that we will refer to as channels. The 10-11 and 
11-12 keV channels have a lower signal-to- noise ratio because of an 
instrumental background line feature. Also displayed in Fig. 1 are the 
GOES 8 observations in the low-energy band, extending the energy 
coverage below the RHESSI range. 

Visual inspection of Fig. 1 yields many events per orbit. There is no 
time of stationary flux. During the total selected observing time of 373 
minutes, we found 64 events, thus 10.3 per hour on average in one solar 
hemisphere. Lin et a/.(1984) reported 10 h~^, Lin et a/.(2001) 5.5 h~^, 
and Shimizu (1995) 26 h~^, all of them observing during higher levels 
of solar activity. 

Some general properties follow immediately from Fig.l: 

— The largest variability occurs in the channels between 5-9 keV. 

— Some events in the 6-9 keV channels are not visible in the GOES 
light curve, being dominated by photons < 3 keV (see e.g. event 
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Figure 1. Top: Light curve observed by GOES 8 in the 1-8 A band (1.6-12.4 keV) of 
soft X-rays. Bottom: RHESSI light curves in one keV channels from 3-4 keV (top) 
to 14-15 keV (bottom). Each channel is multiplied such that it does not overlap 
with others. The time resolution is 20 s. a) Sun light starts at 19:08 UT. h) Sun 
light starts at 03:11 UT and ends at 04:10 UT. c) Sun light starts at 23:40 UT. d) 
Sun light starts at 01:40 UT and ends at 02:40 UT. 



after B in Fig. Id. GOES shows no event larger than AO.l above 
background) . 

— The emission above about 9 keV peaks generally before the emis- 
sion at lower energies. 

— The event duration decreases continuously from 1.6 keV (GOES) 
to 15 keV. 

— In most events there is no emission perceivable beyond about 
12 keV. 
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Figure 2. Spectrum of event A and background B on 2-May-2002 (Fig. Id) observed 
with RHESSI. The event is fitted with a thermal contribution (dotted) and a 
non-thermal component (dashed) having a break point at 8 keV. The background 
is fitted with only a thermal contribution. 



2.2. Photon Energy Spectra 

Two spectra arc shown in Fig. 2. The data were cahbrated using the 
full detector response matrix. The first spectrum was integrated during 
event A shown in Fig.ld from 01:46:40 to 01:47:40 UT, the maximum 
phase at the higher energies. The total background, taken at the time 
of minimum flux before the event, was subtracted. Only the front 
detectors 1, 3, 5, and 9 have been used, having the lowest dropout 
rates in the time interval. Data dropouts are suspected to be caused by 
heavy cosmic rays interfering with the RHESSI electronics. Although 
the software corrects for the deadtime due to these data gaps, they 
introduce additional noise. 

At low energies the spectrum fits well with an isothermal plasma at 
a temperature of 12.1(ib 0.4) MK and an emission measure of 7(ib2) x 
10^^ cm""^. The best fit is shown in Fig. 2, where all normalized residuals 
are below 0.4. The errors originate mainly from the choices of the back- 
ground and the low-energy break point for non-thermal emission. As 
the power-law distribution of the non-thermal component must turn 
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over at low energies, it was approximated by a constant below some 
energy. The quality of the fit is not sensitive the break energy; best fits 
are between 6 and 9 keV. In the case of Fig. 2 different background 
choices changed temperatures in the range of ±0.8 MK, and different 
choices for the low-energy break for non-thermal emission produced a 
range of ±0.7 MK. The standard error was estimated by taking a third 
of the rms sum of the two. Note that the error in temperature and 
emission measure are anti-correlated. 

At the time of writing, the energies below 5 keV are not yet fully 
calibrated. This can be noted in the fact that single collimators give 
a range of temperature values of ±1.7 MK if using only the lowest 3 
keV. The scatter decreases by a factor of 3 for taking a larger range. At 
6—7 keV, and possibly also at 7 8 kcV, the flux is often enhanced. This 
cannot be reduced substantially by the choice of the break point. The 
enhancement could be caused by the 6.7 keV iron line emissions not 
included in the fit. Therefore we did not use 6-8 keV for fitting. Instead 
the break point was put at 8 keV in accordance with other information 
(Figs. 4 and 5). 

In the energy range of 10 - 15 keV, the isothermal model does 
not fit the observations of microflare A (Fig. 2). A power-law photon 
distribution was therefore added. Its fitted slope depends on the break 
point and has an exponent in the range -4.5 ^ 7 ^ -5.7. For the 
choice of the break point at 8 keV, the best fit yields -4.5. The spectrum 
disappears in the background beyond 15 keV. 

Instead of a power-law, the 10 - 15 keV spectrum could also be fitted 
by a second thermal component with higher temperature and smaller 
emission measure. The best fit is achieved for a temperature of 24.7 
MK and an emission measure of 0.14 x 10^^ cm^'^. The contribution of 
the hot kernel would not be discernible below 8 keV and fit all other 
constraints. 

The thermal energy content of an isothermal plasma is 

Eth - sVMVkBT (1) 

It is assumed in Eq.(l) that the electron and ion densities are about 
equal and constant in the volume V. A4 is the emission measure, and 
ks is the Boltzmann constant. Putting in the observed values, we find 
for the thermal content of the microflare (index A) and of the hot kernel 
(index k) 

4.2 X 10l^^A^ [erg] , (2) 
Ei^l.2xW^*VVk [erg] , (3) 
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Thus the hot kernel would contain less energy than the main component 
for Vk<VA. 

The spectrum at the time of the lowest count rate in the orbit 
(marked B in Fig. Id) is also displayed in Figure 2. It is the time interval 
from 02:14:40 - 02:19:56 UT, when no microflare was observed and the 
most recent one had occurred more than 4 minutes before (Fig. Id). For 
this case, the background interpolated from the satellite's night time 
before and after the observations was subtracted. An isothermal model 
having a temperature of 6.4 MK fits well. No non-thermal or hot kernel 
component is above the background level. The emission measure of the 
hemisphere is measured 45 x lO^^cm"^, and its thermal energy content 
is 

5.6x10^471^ [erg] . (4) 

The fact that Eqs.(2) and (5) have similar factors does not mean that 
K, since probably Vb 3> Va- 

Microflare A can be imaged with RHESSI collimators 3 to 9. The 
microflare took place in NOAA Region 9932 (see image and movie 
on CD) and had a nominal FWHM size of 8"xll". It is close to the 
resolution limit, thus Va ^ 3 x lO^^cm^. Eq.(2) yields an upper limit 
on the thermal energy content of 7.2 x 10^ ''erg. A lower limit for the 
electron density of 4.8 x 10^ cm^'^ is derived from emission measure and 
volume assuming a filling factor of unity. Microflare imaging is studied 
in detail by Krucker et a/. (2002). 

The relation between the temperature of the thermal component 
and the emission measure is further explored in Fig. 3. Temperature 
and emission measure have been determined identically to Fig. 2. Using 
the same background and break point that produce similar systematic 
errors, the remaining relative error between measurements is estimated 
about ±0.2 MK, as determined from different fits. The temperature 
maximum occurs between peak fluxes in 12-15 keV and 3-7 keV. The 
product Ta/^M is proportional to the energy content if the volume 
remains constant. The product increases initially until about 01:48:30 s, 
when the 12-15 keV count rate has dropped to the half-maximum level, 
and then remains approximately constant throughout the microflare. 

2.3. Low AND High Energy Photons 

The light curves integrated from 3-7 keV (low energy) and 12-15 
keV (high energy) are shown in Fig. 3. The low-energy band peaks 
at 01:48:10(ib20) UT. The high-energy count rate reaches the prefiare 
background level of 1.8 counts s"! at 01:48:40(±20) UT. Thus the two 
light curves exhibit the Neupert effect (Neupert, 1968) known in regular 
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Figure 3. Top: Temperature (solid) in units of 10 MK, emission measure (dashed) 
in 10~*®cm~^ and the product T\fM (dotted) of event A in Fig. Id (thermal com- 
ponent). Bottom: Light curve of event A in the energy bands 3 8 keV and 11 15 
keV (multiplied by factor of 10). The integration time is 8 s, the data are calibrated 
by the diagonal matrix elements. 
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Figure 4- Spectrogram observed by RHESSI of event A in Fig. Id. The energy chan- 
nels (1 keV) and time bins (8 s) are presented by pixels. A logarithmic scale is 
applied, bright meaning enhanced photon flux. 



flares. It is generally interpreted there as the signatures of different 
emissions relating to a thermal and a non-thermal electron population 
in the scenario of electron-beam driven evaporation. 

The high-energy band peaks 80 seconds before the low-energy band. 
The FWHP duration at 12-15 keV is only 45 s, whereas it is 148 s at 
3-5 keV. 

Further differences between low and high-energy photons appear 
in the spectrogram (Fig. 4). The data are calibrated using a diagonal 
response matrix. The high-energy peak at 01:46:50 UT is clearly visible. 
It takes a straight vertical shape emerging at about 8 keV out of the 
thermal part. The latter has a triangular shape with an apex at about 
01:47:25 UT and 12 keV. 

The timing of photon fluxes at different energies was studied also 
by cross-correlation. The light curves at 4 s resolution were each cross- 
correlated relative to the most variable channel (7-8 keV). The derived 
cross-correlation coefficients revealed a local maximum at zero lag even 
at energies >20 keV, where no solar photons were detected above back- 
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Figure 5. Correlogram of Fig. 4 showing the cross-correlation of relative to the 7-8 
keV channel. Bright pixels represent high correlation coefficients. 

ground. As the zero-lag cross-correlation peak also exists at satellite 
night, it must be an instrumental effect. Most likely it is due to data 
dropouts, except for the auto-correlation in the 7-8 keV channel. The 
effect was eliminated by interpolation between the neighboring lags and 
integration to 8s. 

Figure 5 depicts the result of the cross-correlation in spectrogram 
form (correlogram). Channels with large background, such as 10-11 and 
11-12 keV, have higher noise and yield lower cross-correlation coeffi- 
cients. Again, the difference between low and high energy components 
appears. Most relevant is the change in drift rate of the cross-correlation 
peak. At energies higher than 9 keV, the peak does not drift but forms 
a vertical structure at a constant lag of -25s. Below 9 keV, the peak 
drifts to larger lags (from upper left to lower right), reaching about 
+25s at 3-4 keV. 

2.4. Low-energy Component 

The 25 largest microflares in the sample have been analyzed similarly 
to the above analysis. Their average peak flux in the 7-8 keV band 
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Figure 6. Drift rate of the light curves of the 1 keV channels between 3 and 9 keV 
vs. peak photon flux in the 7-8 keV channel. 



above the preceding minimum is 0.7 photons s~^ cm~^. The average 
FWHM duration is 131 (±103) seconds (standard deviation), and the 
median value is 114 s. We found no correlation between peak flux and 
duration in agreement with Shimizu's (1995) measurements at lower 
energies. 

A relation between peak flux and absolute drift rate of the cross- 
correlation peak is apparent from Fig. 6. The drift rate was determined 
from the cross-correlation coefficients of the 1-keV channels in the 3- 
9 keV range by linear (Fig. 5) of the peak values in energy and time. 
The peaks in each channel are determined from a boxcar smoothing 
of 40 s and spline interpolation. Spline and regression allow a higher 
resolution and the determination of larger drift rates than from the 
pixels in the spectrogram. Nevertheless, the error increases enormously 
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at drift rates beyond about 2 keV/s. Large microflares have slower 
absolute drift rate, which is usually negative (peak drifting from high 
to low energies). Some smaller microflares have large absolute drift. In 
one case out of 25 it was found positive, but the error range includes 
negative values. The anti-correlation between photon flux and drift rate 
is more prominent for events during the same orbit. 

The spectrum of X-ray microflares at energies in the 3-9 keV range 
is well modeled by an isothermal plasma as shown in the example of 
the previous section. The best fitting temperatures in the 25 largest 
and well isolated events are in the range 10.2 - 14.7 MK, with an aver- 
age value of 12.0(±1.2) MK. This is considerably more than reported 
previously (Shimizu, 1995; Feldman et al., 1996). The difference seems 
to originate from Yohkoh/SXT and GOES favoring the 4-7 MK range, 
but RHESSI observing the most energetic thermal photons and thus 
emphasizing the highest temperatures. This is supported by GOES 
observations of the microflare of Fig. 2 (GOES class A5) where the two 
GOES channels having lower energy bounds of 1.6 and 3.2 keV indicate 
a temperature of only 7±1 MK (background subtracted). 

2.5. High-energy Component 

The second component is best observable in 11-14 keV photons. As 
a test of its presence, the <8 keV spectrum was fitted by a thermal 
model, and the 12-13 keV channel compared to the fit. In 6 of the 25 
largest microflares this channel was clearly above the thermal, 2 more 
cases were doubtful. In 4 cases, all at relatively high temperature, the 
12-13 keV flux was within the error bars consistent with the thermal 
model. In the rest, 13 cases, the 12-13 keV flux was not significantly 
above background and could not be classified. 



3. Discussion 

The existence of a second, non-thermal component at high energies 
in microflares is not obvious. We have searched for it in the following 

ways: 

1. The spectrum (Fig. 2) of microflares below 10 keV is dominated 
by a thermal component with a temperature of the order of 10 MK. In 
24% of the cases an enhanced high-energy tail is observed extending 
up to 15 keV. It could also be interpreted by a hot thermal kernel. 

2. The spectrogram (Fig.4) shows more clearly the difference be- 
tween the two components: (i) The low-energy thermal component 
having a triangular shape with a long decay at the lowest photon 
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energies and (ii) a short vertical structure extending to higher energies 
in the rise phase of the former component. 

3. The drift of the cross-correlation peak at low energies is a combi- 
nation of two effects: (i) The maximum of the light curve is more and 
more delayed the lower the energy, and (ii) the shape of the microflare 
becomes asymmetric at lower energies having a longer decay time than 
rise time (Fig.3). Cross-correlation is thus also a method to identify the 
two components and to separate them in energy. 

There are two major indications for a "non-thermal" origin of the 
high-energy component: 

— A hot kernel of thermally emitting plasma also fits the 9-14 keV 
spectrum, but would produce a drift in the correlogram during 
cooling that is not observed in the second component. 

— The second component has shorter duration, occurs in the early 
phase of the thermal component and vanishes at the peak flux 
of the thermal component. It resembles the non-thermal compo- 
nent of regular flares in its relation to the first component known 
as the Neupert effect. The second component thus exhibits the 
characteristics of thick target bremsstrahlung of an electron beam. 

Therefore, we interpret the second component of microflares as non- 
thermal emission. 

An important result is the increase of the emission measure through- 
out the event (Fig.3). It indicates that the amount of high-temperature 
material augments (excluding the unrealistic possibility of compres- 
sion). At the same time the temperature decreases. A surprising char- 
acteristic is the flat top of the product T^/Ai from the time of the 
temperature peak early in the flare until the late decay phase. Under 
the assumption of constant volume, the flat top suggests a constant 
thermal energy content 

Thermal conduction could explain such a behavior if radiative en- 
ergy loss can be ignored. The thermal component may start with a high 
temperature and a small mass. Thermal conduction then adds more 
material to the temperature range observed, but decreases the mean 
temperature at the same time. It is consistent with the observation of 
the thermal component, having a later peak and longer duration at low 
energies. 

Conduction can also explain the drift of the cross-correlation (Fig. 5). 
As the temperature decreases, the count rate in the higher energy chan- 
nels diminishes more rapidly than at low energies where the increasing 
emission measure even delays the maximum flux. The smaller the flare 
size, the faster this takes place. Thus the cross-correlation peak drifts 
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the faster, the smaller the flare. Considering tlic previously reported 
relation between flare size and peak flux (see Introduction), this may 
explain the observed inverse relation between drift rate and peak flux 
(Fig.6). 

Conduction models have been proposed to interpret a Neupert-like 

relationship between higher and lower energy photons (reviewed in Den- 
nis and Schwartz, 1989). Although the second component of microflares 
at high energies could be fitted by the thermal emission of a hot kernel, 
it would not contain enough energy to allow for conductive heating of 
the first component. Therefore, conduction may explain the evolution 
of the first component, but cannot explain the relation between the 
first and second component. 

We propose here that the flare duration at 3-9 keV (thermal com- 
ponent) is not controlled by radiation loss, but by thermal conduction 
distributing the energy to a larger mass reservoir at lower temperature. 



4. Conclusions 

The X-ray spectrum from 3-15 keV of solar microflares in active 
regions has been observed with high energy resolution by the RHESSI 
satellite. The thermal component can be identified in all microflare and 
traced up to some 12 keV in the spectrogram. It can be well separated 
from a second component at higher energies, not visible in all events, 
that peaks simultaneously at all energies, has short duration and occurs 
during the rise phase of the first component. The second component 
can be traced in energy down to about 8 keV. The existence of two 
components excludes a previous interpretation (see Introduction) as 
one non-thermal electron population having a cut-off distribution in 
energy. 

The similarities of microflares to regular flares suggest identical 
physical processes and similar plasma conditions. The properties of the 
thermal and non-thermal components are consistent with the standard 
flare scenario of precipitating electrons heating cold material to flare 
temperature. The collisional energy loss of beaming electrons causes 
bremsstrahlung emission observable in the 9-15 keV range and heats 
cold material to emit thermally at <12 keV. 

High spectral resolution at the transition between thermal and non- 
thermal emission has been found to be useful for studying the role 
of energy deposition by beams and thermal conduction. No-shutter 
observations by RHESSI are well suited to further quantify these re- 
sults. This will lead to firm estimates on the various forms of energies 
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involved in microfiares and determine their role in heating the active 
region corona. 
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